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Abstract: Garra imberba is widely distributed in China. At the moment, both Garra yiliangensis and С. hainanensis are treated as 
valid species, but they were initially named as a subspecies of G pingi, a junior synonym of С imberba. Garra alticorpora and С; 
nujiangensis also have similar morphological characters to G imberba, but the taxonomic statuses and phylogenetic relationships of 
these species with С: imberba remains uncertain. In this study, 128 samples from the Jinshajiang, Red, Nanpanjiang, Lancangjiang, 
Nujiang Rivers as well as Hainan Island were measured while 1 mitochondrial gene and 1 nuclear intron of 24 samples were 
sequenced to explore the phylogenetic relationship of these five species. The results showed that С: hainanensis, G. yiliangensis, С 
alticorpora and С. imberba are the same species with С. imberba being the valid species name, while G nujiangensis is a valid 


species in and of itself. 
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With 105 valid species Garra is one of the most 
diverse genera of the Labeoninae, and has a widespread 
distribution ranging from East Asia to Africa (Froese & 
Pauly, 2012). Menon (1964) first revised the genus and 
divided it into 4 groups and 9 complexes. Garra 
imberbis and G. imberba were classified into the 
imberbis complex, which is distinguished from the other 
complexes by having more lateral line scales, a larger 
body size, a shorter distance between vent and the pelvic 
fin base, and no barbels (Menon, 1964). 

Garman (1912) described Garra (Ageneiogarra) 
imberba using samples from the Min River, Kiating, 
Szechuan (Leshan, Sichuan), China. С imberba can be 
distinguished from other Garra species by its lack of 
barbels, 46-52 lateral line scales, 16 circumpeduncle 
scales, no secondary rostrum, and an anus close to pel- 
vic-fin base. Ageneiogarra was treated as a synonym of 
Garra, because the number of barbels was thought to be 
not eligible to define a new subgenus during the initial 
intensive study of cyprinid fish. (Karaman, 1971; 
Kottelat, 1998). Fang (1943) later renamed the materials 
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from Jinshajiang River as С. pingi, but Kottelat (1998) 
treated С: imberba as a valid name, and this convention 
is still widely accepted at present. 

Four subspecies of Garra imberba have also been 
identified, although their species validity is still uncertain. 
Wu et al (1977) described a new subspecies, G pingi 
yiliangensis, based on the materials from Nanpanjiang 
River. G. pingi yiliangensis has 9 branched dorsal-fin 
rays, a standard length 5.7—6.1 times body depth, and the 
caudal-peduncle is 2.0–2.1 times longer than its depth. 
Zheng & Chen (1983) described another subspecies from 
Hainan Island, named С. pingi hainanensis. It has 46—47 
lateral-line scales, and the length of caudal-peduncle is 
1.1-1.2 times its depth. Garra imberba, G. рта 
yiliangensis, and С. pingi hainanensis all have isolated 
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distributions, so they were treated as distinct species 
under the phylogenetic species as conceptualized by 
Kottelat (1998). Garra pingi caught from the Nujiang 
River (upper Salween River) were later described as С 
nujiangensis by Chen et al (2009). It has 12—14 circump- 
eduncle scales, and some of its characteristics are meri- 
stically and metrically different from С imberba. 
However, Chen et al also found that the new species has 
some characteristics similar with С. imberba, especially 
among juveniles, so they just temporarily treated it as a 
distinct species from С. imberba. Lastly, Garra 
alticorpora was named by Chu & Cui (1987) based on 
two specimens from Pingbian, Yunnan, which belongs to 
the Red River Drainage system. The most significant 
characteristic of this Garra is that its body depth is 
longer (vs. shorter) than its head length. After original 
sampling, however, specimens collected from the same 
location did not have the same characteristics, so the two 
G. alticorpora samples were thought to be gravid females 
and treated as a synonym of G imberba by Zhou et al 
(2005). Besides С. imberba, all four of these species (С. 
pingi yiliangensis, С. pingi hainanensis, G. nujiangensis, 
G. alticorpora) were known as either subspecies of it or 
considered synonym of it, and their species validity and 
phyloge netic relationship is still uncertain. 


Hainan Island 
Red River 


Mekong River : 
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The mitochondrial DNA and nuclear gene have been 
widely used in the studies of molecular phylogenetic 
analyses of subfamily Labeoninae or genus Garra (He et 
al, 2008; Li et al, 2005; Tang et al, 2009; Yang & 
Mayden, 2010; Zheng et al, 2010; Yang et al, 2012). In 
this study, morphological and molecular methods were 
used to investigate the phylogenetic relationships among 
the five species or subspecies of G imberba and its 
related species distributed in China. 


MATERIALS AND METHODS 


Sampling 

A total of 128 specimens were measured for 
morphologic analysis, including 90 Garra imberba from 
Jinshajiang and Red-Mekong River, 9 G yiliangensis 
from Nanpanjiang River, 11 С; hainanensis from Hainan 
Island, 16 С. nujiangensis from Nujiang River and 2 С 
alticorpora from Red River. A total of 24 samples were 
used in molecular phylogenetic analyses. The collection 
localities are shown in Figure 1. Garra yiliangensis is 
only distributed in Nanpanjiang and Beipanjiang 
Rivers(upper Pearl River), but it was not included in this 
study for molecular analysis because no specimen has 
been caught in the last decade. No new G alticorpora 














Figure 1 Geographic distribution of samples used in molecular analyses 
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specimens have been caught since it was named, so it 
was also not included in the molecular analysis. In 
addition, one sample of Garra fasciacauda from Pwer, 
Yunnan was sequenced. This species was thought only to 
be distributed in the middle of the Mekong River, such 
as in Thailand and Cambodia, and this is the first record 


Table 1 Taxa used in molecular analyses with source of each 


of its presence in China. 

All the specimens used in this study are deposited in 
the Collection Room of Fishes, Kunming Institute of 
Zoology (KIZ), Chinese Academy of Sciences (CAS). 
Voucher information and GenBank accession numbers 
are listed in Table 1. 


samples and GenBank accession numbers for each gene 


GenBank Accession No. 

















No Drainage Voucher specimen Taxon Location 
cyt b S7 
1 YZ1 KIZ041110093 Garra imberba Chuxiong, Yunnan, China KC119045 KC119071 
2 YZ2 KIZ2010001891 Garra imberba Huize, Yunnan, China KC119046 KC119072 
3 RRI KIZ2010010001 Garra imberba Xinping, Yunnan, China KC119047 KC119073 
4 RR2 KIZ2010010002 Garra imberba Hekou, Yunnan, China KC119048 KC119074 
5 LX1 KIZ2004001193 Garra imberba Puer, Yunnan, China KC119049 KC119075 
6 LX2 К172009002070 Garra imberba Lvchun, Yunnan, China KC119050 KC119076 
7 МК! К172008000088 Garra imberba Xishuangbanna, Yunnan, China KC119051 KC119077 
8 MK2 К172008008269 Garra imberba Puer, Yunnan, China KC119052 KC119078 
9 НМ К172008003973 Garra hainanensis Ledong, Hainan, China KC119053 KC119079 
10 HN2 К172008003938 Garra hainanensis Ledong, Hainan, China KC119054 KC119080 
11 SW1 KIZ2005005142 Garra nujiangensis Cangyuan, Yunnan, China КС119055 КС119081 
12 $\/2 KIZ2005005145 Garra nujiangensis Cangyuan, Yunnan, China KC119056 KC119082 
13 KIZ2008005732 Garra micropulvinus Wenshan, Yunnan, China KC119057 KC119083 
14 KIZ2011002800 Garra findolabium Jinping, Yunnan, China 10864598 KC691274 
15 KIZ5171 Garra cryptonemus Liuku, Yunnan, China JQ864587 KC691275 
16 KIZ2007002789 Garra caudofasciatus Jiangcheng, Yunnan, China 10864588 KC691276 
17 К172004014928 Garra fasciacauda Puer, Yunnan, China JQ864597 KC691277 
18 KIZ2006004460 Garra tengchongensis Tengchong, Yunnan, China 10864586 KC691278 
19 KIZ2004000816 Garra dulongensis Dulongjiang, Yunnan, China 10864590 KC691279 
20 KIZ2006004422 Garra qiaojiensis Yingjiang, Yunnan, China 10864583 KC691280 
21 КІ22006003543 Garra salweenica Dehong, Yunnan, China 10864593 KC691281 
22 KIZ2005000086 Garra orientalis Longlin, Guangxi, China 10864581 KC691282 
23 KIZ2005002392 Garra mirofrontis Yunxian, Yunnan, China JQ864584 KC691283 
24 Labeo senegalensis AB238968 AY 103160 


Morphological analyses 

Measurements were taken point to point with a 
digital caliper to the nearest 0.1 mm. Measurements and 
counts were made on the left side of individuals when- 
ever possible (Kottelat, 2001). Predorsal, prepectoral, 
prepelvic and preanal lengths were taken respectively, 
from the anterior most tip of the snout to the dorsal-, 
pectoral-, pelvic- and anal-fin origins. Interorbital width 
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was measured between the upper margins of the eyes. 
Abbreviations used in this paper are as follows: SL 
(standard length), HL (head length), HD (head depth), 
HW (head width), DFL (dorsal fin length), PDL (pre- 
dorsal fin length), PFL (pectoral fin length), PPL 
(prepectoral fin length), VFL (ventral fin length), 
PVL(preventral fin length), AFL (anal fin length), PAL 
(preanal fin length), CFL (caudal fin length), CPL(caudal 
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peduncle length), CPD (caudal peduncle depth), DAA 
(distance between anal and origin of anal fin), DVA 
(distance between ventral fin base and origin of anal fin), 
DL (disk length), DW (disk width), ED (eye diameter), 
IOW (interorbital width), SNL (snout length). 

Principal component analysis (PCA) was used to 
examine the significance of difference among samples 
of different species or subspecies with SPSS 13.0. 
Default settings of factor analysis were applied for 
data standardization. Covariance matrix was used in 
the analysis. All samples were assigned into 5 groups 
based on species. 


DNA extraction and amplification 

Total genomic DNA was extracted from Фе 
ethanol preserved specimens using а standard 
phenol/chloroform isolation and ethanol precipitation. 
Polymerase chain reaction (PCR) was used to amplify 
the sequences of one mitochondrial gene and one nuc- 
lear gene: cytochrome b (cyt b) and the first intron of 
S7 ribosomal protein (r-protein) gene. These genes are 
considered to be effective for understanding the evolu- 
tionary relationships among populations and species 
(He et al, 2008; Irwin et al, 1991). The complete cyt b 
gene was amplified with primers L14724 and H15915 
(Xiao et al, 2001), S7 was amplified using primers 
57КРЕХІЕ апа S7RPEX2R (Chow and Hazama, 1998). 
The PCRs contained approximately 100 ng of template 
DNA, 1 uL of each primer, 5 uL of 10x reaction buffer, 2 
uL dNTPs (each 2.5 mM), and 1 U Taq DNA polymerase 
in total 50 uL volume. The PCR profile consisted of an 
initial denaturation step (3 min at 94 °C), followed by 35 
cycles performed in the following order of denaturation 
at 94 °С for 1 min; annealing at 52 °С for 1 min; and 
elongation at 72 °С for 1 min; and a final extension at 
72°C for 8 min. PCR amplification products were 
purified using the Sangon DNA purification kit 
according to manufacturer’s instructions. The purified 
PCR product was sequenced by the Shanghai Sangon 
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Sequence alignment and phylogenetic analyses 

DNA sequences were edited using DNASTAR 
v7.1.0 (DNASTAR Inc.). All sequences were aligned by 
Clustal W and checked by eye in MEGA 5.0 (Tamura et 
al, 2011). Amino-acid translation was also conducted 
using MEGA 5.0 (Tamura et al, 2011) to confirm the 
correct reading frame positions and find unexpected stop 
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codons. The genetic distances (p-distance with 1 000 
bootstrap) of the two sequences between taxa were 
calculated using MEGA 5.0 (Tamura et al, 2011). Two 
operational datasets were constructed for subsequent 
analyses: (i) mitochondrial DNA sequences, complete 
cyt-b (1 140 bp); (ii) nuclear DNA sequence, only 
include the first intron of S7 (805-809 bp). 

Phylogenetic analyses were performed using maxi- 
mum likelihood estimation (ML), and Bayesian inference 
methods with PAUP* version 4.0 510 (Swofford, 2003) 
and MrBayes 3.1.2 (Ronquist & Huelsenbeck, 2003) 
respectively. For the ML analyses, a heuristic search was 
adopted with 1000 bootstrap replications, and a random 
tree was used as starting tree. Bayesian analyses was 
conducted for 5 million generations, using one cold and 
3 incrementally heated chains and sampling every 100 
generations. Two independent runs were performed to 
confirm convergence. Models were calculated by jMod- 
eltest v0.1.1 (Posada, 2008). 


RESULTS 


Morphological analyses 

Counts and proportional measurements of Garra 
imberba, С. yiliangensis, G. hainanensis, С. nujiangensis 
and G. alticorpora are shown in Table 2. The factor 
analysis results showed that the first 3 components can 
explain 73.8% of the total variances. For the first 3 
components, a scatter graph was drawn as follows (Figure 
2). The 5 biggest correlations of the factor analyses of the 
first 3 components is shown in Table 3. Factor 1 reflects a 
primarily more slender body shape, and is composed 
mostly of samples of G. yiliangensis that are not clustered 
with the other studied species. There was no significant 
difference for factor 2 and 3 among these species. 


Molecular analyses 

Among the two sequences used in this study, only 
cyt b is a protein coding gene and could be translated 
into amino acid sequences without interruption. The 
pairwise distances of cyt b and S7 are shown in Tables 
4, 5 and 6. The overall mean p distance of clade Al (G 
imberba and С. hainanensis ) based on cyt b and S7 are 
1.9%-0.3% and 0.2%+0.1%. species of clade АЗ 
contains valid species of Garra distributed in China and 








the mean distance within this clade based on cyt b and 57 
аге 12.3%+0.5% and 11.2%+0.7%. The significant 
difference within clade Al and A3 implied that samples 








of clade A1 actually represent the same species. 
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Figure 2 Scatter plots analyzed by PCA for samples from different locality groups based on 22 morphometric characters 

A: scatter plot of scores оп factorl and factor 2; В: scatter plot of scores on factor2 and factor3; AL: С alticorpora, from Red River; 
HN: С hainanensis, from Hainan Island; LX: С imberba, from Lixianjiang, the biggest tributary of Red river; MK: G. imberba, from 
Lancangjiang River; NP: G yiliangensis, from Nanpanjiang River, the upper stream of Pear River; RR: G imberba, from the main 
stream of Red River; SW: G nujiangensis, from Nujiang River; YZ: G imberba, from Jinshajiang River. 


Table 3 Proportion characteristics and principal components extracted from examined materials 





No. Component 1 Correlation Component 2 Correlation Component 3 Correlation 
1 0.966 CPD/CPL 0.842 DW/HL 0.648 IOW/HL 
2 —0.859 CPL/SL 0.773 HW/HL —0.457 DL/HL 
3 0.648 HL/SL 0.577 DL/HL —0.45 PPL/SL 
4 0.606 HD/SL 0.408 IOW/HL —0.405 HL/SL 
5 0.602 PPL/SL 0.352 CPD/SL —0.404 DAA/DVA 
The BI and ML trees of cyt b and S7 have the same 

topology. Garra imberba from Jinshajiang, Red, Mekong DISCUSSION 


River and G hainanensis are clustered together, and 
Garra nujiangensis is the sister species of С: imberba. In 
general, phylogenetic relationships based on the phylog- 
enetic trees (Figure 3) of mtDNA exhibit a marked 
pattern based on geographical distribution, which is 
different from the result of morphological analyses 
(Figure 2). Garra hainanensis has a closer relationship 
with G. imberba from Jinshajiang and the main stream of 
Red River, and then clustered with G imberba from 
Lixianjiang and Lancangjiang River materials. The 
results show a conflict in the modern river systems; 
Lixianjiang River is the biggest tributary of Red River, 
but samples from Lixianjiang have a closer relationship 
with Lancangjiang River. 

G. nujiangensis from Nujiang River can also be 
easily identified (Figure 4) based on the phylogenetic 
trees based on the S7, and a mix of samples from 
Jinshajiang, Red, Lixianjiang, Lancangjiang River and 
Hainan Island. 
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Morphological discrimination between species 

Morphological characteristics are still the most 
important evidence for defining a new species or subs- 
pecies in cyprinid fish. Such analyses usually include 
two kinds of data: meristic and metric. A difference in 
meristic characteristics usually can be used as evidence 
for defining a new species; metric data usually is used to 
describe a new species or subspecies based on significant 
difference (Chu & Cui, 1987; Chu & Chen, 1989). 

Table 2 shows the details of the morphological 
characters of these five species. Most of the meristic data 
do not have significant difference, except for the 
circumpeduncle scales: С. imberba usually has 16 or on 
rare occasions 18 (1 in 90) circumpeduncle scales; most 
of С. nujiangensis has 12, and only 1 sample has 14 
circumpeduncle scales. Number of circumpeduncle 
scales was treated as very important evidence in defining 
species, such as С: giaojiensis and G. orientalis (12 vs. 16) 
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Table 6 Overall mean p distance within and between clades 





computed by MEGA 
Sequence Clade Mean distance (%) SE (%) 
Al 1.9 0.3 
A2 0 0 
A3 123 0.5 
cyt b 

Al Vs. A2 9.5 0.8 

А1 Уз. АЗ 12.7 0.6 

A2 Vs. A3 12.4 0.6 
Al 0.2 0.1 
A2 0 0 
A3 11.2 0.7 

S7 

Al Vs. A2 1.4 0.4 

Al Vs. A3 9.7 0.7 

А2 Уз. АЗ 9.6 0.7 

(Chu & Cui, 1987). Garra nujiangensis сап ђе 


distinguished from all the other species by having 12 
circumpeduncle scales (vs. 16) which implies it is a valid 
species and distinct from G imberba. 


1/100 


з 


The number of branched dorsal fin rays is thought to 
be a stable characteristic and efficient in defining new 
species (Chu & Cui, 1987). In this study, we found it is 
actually not as stable as the circumpeduncle scales. For 
example, in С. imberba, among the 90 examined 
materials 2 samples have 10 branched dorsal fin rays, 10 
have 8, and 72 have 9. The number of lateral line scales 
usually has a changing range and was treated as a general 
characteristic in defining species (Chu & Cui, 1987). 
Meristic data in Table 2 show that these two characte- 
ristics have little but not significant difference among G 
imberba, G. yiliangensis and G. hainanensis, so it is not 
reasonable to treat them as distinct species based solely 
on the above characteristics. 

In the PCA analysis, distinction among these five 
species is not very clear. From Figure 2, it can be seen 
that G. yiliangensis has the most significantly different 
from the other species, but still has 1 material mixed 
within it. Component | mainly stands for the CPD/CPL, 
which means most С. yiliangensis have a slender body, 
but the metric character is not reliable in defining species 
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Figure 3 BI tree constructed using MrBayes based on cyt b gene 
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Figure 4 BI tree constructed using MrBayes based on the 57 


especially in genus Garra (Chu & Cui, 1987). The 
morphological characteristic is closely related with the 
environment, for example eyes will degenerate if the fish 
live in dark caves for a long time. Garra yiliangensis has 
a slender body especially for the caudal peduncle. From 
an ecological view, it is a characteristic adopted for fast 
flowing waters. But this difference is not significant, 
which implies that G. yiliangensis from Nanpanjiang 
River only recently became isolated from adjacent rivers, 
so it should have a close relationship with G imberba 
from Jinshajiang and Red River. 

Garra alticorpora can only be distinguished from 
other related species by the fact that its body depth is 
longer than its head length. There is no difference in 
meristic characteristics with С. imberba and mixed 
with the others in PCA analysis. The only one 
holotype and one paratype were collected in Pingbian, 
Yunnan, which belongs to the Red River Drainage 
system. In the following surveys, however, only С. 
imberba were caught, so the two specimens of G 
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alticorpora may be just a special case or gravid 
females and should be treated as a synonym of G 
imberba (Zhou et al, 2005). 

From traditional morphological characteristics 
analyses, only С. nujiangensis can be distinguished 
from other species. Garra imberba, G. yiliangensis, G. 
and С 
indifferent from one another, and should be assigned 


hainanensis alticorpora ате relatively 
to the same species. Because of the geographical 
isolation, G. yiliangensis and G. hainanensis should be 
treated as different geographical populations of G 


imberba. 


Molecular phylogenetic relationship 

It is widely accepted that Jinshajiang, Nanpanjiang 
and Red River have a close relationship; indeed, 
Lancangjiang, Nujiang and Irrawaddy River have a 
close relationship (Chu, 1986). Based on geological 
studies, sea-level changes during the Pleistocene may 
have reached the maximum lows, which means the 
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Red River might have been linked together with rivers 
on Hainan Island (Hanebuth et al, 2000; Voris, 2000). 
Therefore, С; hainanensis and С. imberba from Red and 
Jinshajiang River possibly have the same ancestor, and 
split only after the Pleistocene era, which may explain 
why С. hainanensis has a closer relationship with С. 
imberba from Red and Jinshajiang River, and why 
they were clustered in a lineage of Gl. The close 
relationship coincides with the geographical process 
of the three rivers. 

The Mekong River has the highest fish diversity 
in Southeast Asia and its fish fauna is also different 
from the Red, Jinshajiang and Nanpanjiang Rivers 
(Froese & Pauly, 2012; Yap, 2002; Zakaria-Ismail, 
1994). Lixianjiang (Song Da in Vietnam) is the 
biggest tributary of the Red River, which originated 
from the Wuliang Mountain, Yunnan, China and flows 
into the Red River at Viet Tri, Vietnam, but the 
mtDNA 
Lixianjiang River have a closer relationship with 
Lancangjiang than the Red River. Zhou et al (2010) 
compared the morphometrics of Vanmanenia tetraloba 


analysis shows that, specimens from 


from Lancangjiang, Lixianjiang and main stream of 
Red River. Their results show that specimens from 
different rivers actually represented different species. 
Garra imberba is not a fish that migrates long 
distances, and so a gene exchange between the main 
stream of Red and Lixianjiang would be very difficult 
given the distance. Therefore, С. imberba from 
Lixianjiang and the main stream of the Red River 
would not have exchanged genes for a long time, and 
they were not clustered in the same clade. Clade Al 
can be divided into two lineages, G1 and G2, which 
coincide with the geographical process and the fish 
fauna of the rivers but conflict with the morphological 
classification of С. imberba and related species. 
Furthermore, the mean p-distance within the clade Al 
is 1.9%, and the biggest p-distance within clade A1 is 
3.2%, which is between С hainanensis and С 
imberba from Lancangjiang River. This value is 
smaller than the minimum value within the genus 
Garra (Yang et al, 2012). With respect to the 
evolutionary history and phylogenetic relationship of 
the samples from different rivers, it is better to treat 
them as the same species. 

The Ist intro of S7 ribosomal protein gene was 
thought to be an efficient marker for the phylogenetic 
analyses of subfamily level (He et al, 2008). It has a 


Zoological Research 


faster evolution rate than the coding nuclear gene but a 
slower rate than mtDNA, so we used it here to explore 
the relationship at a species level. Based on the NJ tree 
of S7 data, G nujiangensis from Nujiang River and G 
micropulvinus can be easily identified (Figure 4), which 
implies С. nujiangensis is a valid species and 57 gene is 
an efficient means for exploring phylogenetic relation- 
ships within genus Garra. Garra imberba and G. haina- 
nensis were mixed together and the mean p-distance 
within this clade is 0.2%. It is obvious that G imberba 
and а: hainanensis are the same species and distinct from 
G. nujiangensis. 


CONCLUSION 


For the morphological analysis, a stable meristic 
characteristic can be treated as an important proof when 
defining a new species, such as circumpeduncle scales, 
but metric characteristics should be considered carefully, 
especially in exploring relationships between species. 
Molecular phylogenetic analysis is an efficient method in 
estimating the phylogenetic relationships between diffe- 
rent species or geographical groups. To make a reason- 
able taxonomy, the morphological differences, phyloge- 
netic relationship and the biogeographic processes sho- 
uld all be considered. 

In this study, both morphologic and molecular 
analyses testified that С; nujiangensis is a valid species 
and distinct from G imberba. Samples from Jinshajiang, 
Red, Lixianjiang, Nanpanjiang, Lancangjiang River and 
Hainan Island have the same ancestor, which was 
isolated because of the formation of modern river 
Different 
morphological 


systems. environments 
different 


differences between these materials is not enough to 


usually lead to 
characteristics, but the 


define new species. Based on the morphological and 
molecular analyses, the samples from Jinshajiang, 
Nanpanjiang, Red, Lancangjiang River and Hainan 
Island have a close phylogenetic relationship and 
should be treated as the same species, and because of 
the geographical isolation, С vyiliangensis and С 
hainanensis should be treated as different geographical 
populations of G imberba. 
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Materials examined 
Garra imberba 

Lancangjiang River: KIZ200800826-829, 3, 
97.8-182.1mm SL, Zhaji River, Jingfu Xiang, Jingdong, 
Yunnan; KIZ2002008504-8512, 9, 60.5-103.3 mm SL, 
Mengga River, Jinggu, Yunnan; KIZ2008000087-89, 2, 
83.1—128.6 mm SL, Menggen county, Lancang, Yunnan; 
KIZ2008008269, 1, 97.6 mm SL, Puer, Yunnan. 

Main stream of Red River: KIZ1964000024-25, 2, 
14.7–16.2 mm SL, Yuanjiang, Yunnan; KIZ2007002528, 
1, 19.4mm SL, Hekou, Yunnan; KIZ2008000922, 
KIZ2008000927, KIZ2008000929, KIZ2008000932, 
K1IZ2008000936, KIZ2008000938-940, 8, 15.1—20.7 mm 
SL, Namoguo, Wenshan, Yunnan; KIZ2010010001, 1, 
102.9 mm SL, Xinping, Yunnan; KIZ2010010002, 1, 
71.6 mm SL, Simao, Yunnan. 

Jinshajiang River: KIZ1982001210-1224, 15, 
47.4—118.85 mm SL, Yanjin, Yunnan; KIZ2004011252- 
253, 2, 69.8—195.0 mm SL, Yongsheng, Lijiang, Yunnan; 
K1IZ2008006836, KIZ2010001891, 2, 101.5—136.33 mm 
SL, Huize, Yunnan; KIZ041110093-99, 6, 91.2—150.7 mm 
SL, Yongren, Yunnan; KIZ1977000732-734, KIZ1977- 
000737-738, 5, 51.1-73.1 mm SL, Недіпе, Yunnan; 
KIZ2008006837-40, 4, 108.5-159.3 mm SL, Huize, 
Yunnan. 

Lixianjiang River: KIZ2000001439-1454, 16, 
64.5—163.1 mm SL, Jinping, Yunnan; KIZ2007002903- 
2906, 4, 68.3—107.2 mm SL, Jiangcheng, Yunnan; 
KIZ2008002405-2416, 12, 102.0-130.6 mm SL, Simao, 
Yunnan; KIZ2004001193, 1, 82.0 mm SL, Puer, Yunnan; 
KIZ2009002070, 1, 53.8 mm SL, Lvchun, Yunnan. 
Garra yiliangensis 

KIZ1960000569, KIZ1963000385, KIZ197700093- 
3-934, 4, 146.1-187.0 mm SL, Yiliang, Yunnan; 
KIZ1977000954-958, 5, 166.3–252.5 mm SL, Luoping, 
Yunnan. 

Garra hainanensis 
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